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We present a numerical solution of the four-quark Bethe-Salpeter equation for ground-state scalar 
tetraquarks with = 0'*’’''. We find that the four-body equation dynamically generates pseudoscalar- 
meson poles in the Bethe-Salpeter amplitude. The resulting tetraquarks are genuine four-quark states that 
are dominated by pseudoscalar meson-meson correlations. Diquark-antidiquark contributions are sub¬ 
leading because of their larger mass scale. In the light quark sector, the sensitivity of the tetraquark wave 
function to the pion poles leads to an isoscalar tetraquark mass Ma ~ 350 MeV which is comparable to 
that of the cr//o(500). The masses of its multiplet partners k and ao/fo follow a similar pattern. This pro¬ 
vides support for a tetraquark interpretation of the light scalar meson nonet in terms of ’meson molecules’. 


PACS numbers: 14.40.Be, 14.40.Rt, Il.lO.St, 12.38.Lg 

Introduction. — The nature of the light scalar mesons 
(t//o(500), k(800) and Qq, /o(980) has been an ongoing 
conundmm for several decades. Although the u is a cen¬ 
tral ingredient to hadronic interaction models, even its ex¬ 
istence was experimentally not firmly established until re¬ 
cently. Taking into account evidence from new data for 
heavy-meson [1] and kaon decays [2], together with disper¬ 
sive approaches to tttt scattering using Roy equations [3], 
the <7 is now again listed in the PDG with a T— matrix pole 
position (400 ... 550) — i(200 ... 350) MeV [4]. 

Nevertheless, the cr and its multiplet partners still do 
not fit well into the light qq meson spectrum. In the non- 
relativistic quark-model classification the pseudoscalar and 
vector mesons = 0 '',1 )ares waves without or¬ 
bital angular momentum, whereas the O^'*", I"*"”, and 
2'^^ states are p waves and therefore their masses should 
be considerably larger, in contrast to what is found exper¬ 
imentally for the scalars. Even more puzzling is the mass 
degeneracy of Oq and /q inside the multiplet. If they were 
ideally mixed qq states, the isosinglet a and isotriplet Oq 
would be mass-degenerate and consist of light quarks only, 
whereas k would contain one strange quark and /q two. 
A qq description would then require an enormous amount 
of flavor mixing to generate the observed splitting between 
a and Og and, in addition, to establish an accidental mass 
degeneracy between Og and /g. Connected with this are 
the different decay properties: Og and /g are relatively nar¬ 
row states close to the KK threshold, whereas a and k are 
broad resonances that decay into vrvr and vriT, respectively. 

A natural explanation for these peculiarities was pro¬ 
posed long ago [5]: what if the light scalar mesons were 
tetraquarks in the form of diquark-antidiquark (dq-dq) 
states? A scalar diquark forms a color-5f7 (3)c and flavor¬ 
s' f7 (3)/ anti triplet, and the combination 3 (8> 3 = 1 © 8 
provides a color singlet and a flavor nonet. However, 
the mass ordering is reversed: the a would be the light¬ 
est state made of u/d quarks only, whereas /g and Og 
would be heaviest and mass-degenerate because they carry 
two strange quarks. This would also explain the decay 
widths: /g and Ug are close to KK threshold and there¬ 


fore narrow, but a and k can simply fall apart through 
the gluonless OZI-superallowed decay. In such a scenario 
the true scalar qq ground states could be identified with 
the experimental ‘first excited nonet’ with masses in the 
1.3 ... 1.5 GeV region, similar to the axialvector and ten¬ 
sor mesons and in agreement with the nonelativistic esti¬ 
mate. The non-qq interpretation of the light scalar nonet 
is supported by a variety of theoretical approaches such as 
QCD sum rules [10], unitarized ChPT [11, 12], relativis¬ 
tic quark models [13], effective Lagrangians [14], or the 
extended linear a model [15]. 

Of course, the four quarks may arrange themselves also 
differently: apart from the dq-dq configuration, they could 
form meson molecules with two color-singlet clusters. In 
the heavy quark region, tetraquark candidates in the XYZ 
spectrum such as the X(3872) have been suggested to be 
of such molecular nature [8, 9]. In the light meson sec¬ 
tor, the decay patterns and the proximity of Og and /g to 
the KK threshold may also point in this direction [6, 7]: 
the KK molecule picture naturally explains the observed 
narrowness of these states. A potential drawback of the 
molecule picture for the light meson sector is that in poten¬ 
tial models for light scalar molecules only some states are 
found to be bound, so that no complete nonet exists [6]. 

In this work we reconsider the question of the tetraquark 
nature of scalar mesons. We solve the relativistic four-body 
Bethe-Salpeter equation (BSE) for two quarks and two an¬ 
tiquarks with the overall quantum numbers of a scalar. 
Compared to previous approaches to the four-body system 
we improve on the following aspects: (a) An inherent limi¬ 
tation of many dq-dq and molecular models is that they as¬ 
sume specific quark configurations inside the tetraquark to 
make the model tractable. Depending on the basic assump¬ 
tions on the dominating part of the underlying QCD forces 
these are organized either in terms of dq-dq or molecular 
flavor states, and corresponding assumptions on the spatial 
part of the wave function are made. We overcome such lim¬ 
itations by solving the four-quark BSE without any preju¬ 
dice on the flavor and spatial structure of the wave func¬ 
tion. (b) The framework is quantum held theoretical in 
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nature and fully relativistic in contrast to potential mod¬ 
els. (c) The (anti-)quarks acquire a dynamical, momentum- 
dependent mass via their nonperturhative interactions as 
described hy the Dyson-Schwinger equation of the quark 
propagator. The quark-(anti-)quark interaction in our ap¬ 
proach is given hy nonperturhative gluon exchange in an 
approximation that satisfies the axial Ward-Takashi iden¬ 
tity. As an important consequence, the Goldstone nature of 
the pseudoscalar mesons is preserved. 

It turns out that especially the last property is crucial 
for a successful description of the light scalar meson spec¬ 
trum. Driven hy the underlying quark-gluon interaction, 
meson and diquark pole structures are generated dynami¬ 
cally hy the two-hody interactions in the four-hody BSE. 
Due to the color algebra, the interaction in the qq channel 
is by a factor of two stronger than the one in the qq chan¬ 
nel. Consequently, diquarks have a larger mass scale than 
corresponding mesons with the same flavor content. In 
the light quark sector, dynamical chiral symmetry breaking 
greatly enhances this difference. As will be explained in 
detail below, these light meson poles dominate the resulting 
tetraquark wave function whereas the diquark singularities 
only play a very minor role. Thus a ’molecular picture’ of 
light scalar mesons naturally arises. We obtain a complete 
multiplet of states including the /o(980) and ao(980) with 
strong iC.^-components, but also the /o(500) with a dom¬ 
inating TTTT component which naturally explains its large 
decay width. 

While we arrived at part of these conclusions already in 
the two-body framework of Ref. [16], the current work is 
a substantial improvement since it does not rely on many 
of the assumptions made in [16]. Still, there are approxi¬ 
mations involved to make the extremely complicated equa¬ 
tions tractable. These are explained and discussed below 
in the technical part of the paper. We will argue that all 
approximations affect the quantitative but not the qualita¬ 
tive results of this work. In particular we believe that the 
molecular interpretation of the states is a robust feature of 
the framework. 

In order to avoid confusion there is an important caveat 
to make: Here and in the following we use the term 
’molecule’ or ’molecular picture’ in a probably different 
sense than what is traditional in the literature. Working 
in momentum space, we do not make any assumptions 
on the spatial distribution of the (anti-)quarks inside the 
tetraquark. Thus we cannot, and in fact do not aim to, 
make a distinction between a state with spatially sepa¬ 
rated (anti-)diquark or meson clusters. What distinguishes 
a ‘diquark-antidiquark’ from a ‘meson molecule’ picture in 
our framework is solely the influence of the corresponding 
analytic structure in the internal bound-state kinematics on 
the tetraquark wave function (as detailed below). 

While we exemplify our approach with the important 
case of scalar tetraquarks, the framework is general and 
has potential applications also in many other branches of 
theoretical physics. Obvious examples are four-body states 




FIG. 1. Four-body BSE for the tetraquark BS amplitude. Only 
two-body interactions are retained, whereas irreducible three- 
and four-body interactions are neglected. 


on the nuclear or the atomic level with interesting (non- 
relativistic) applications pointed out recently [17]. 

This work is organized as follows: In the next two sec¬ 
tions we outline the theoretical framework and detail the 
approximations that we introduce in order to make the four- 
body equation tractable. We then present our results and 
discuss some implications in the conclusions. 

Bethe-Salpeter equation. — If four-quark states exist 
in QCD, they will appear as poles in the qqqq scattering 
matrix. The T matrix satisfies an exact scattering equation 

T = K + KGoT, (1) 

where K is the four-quark interaction kernel and Gq is the 
product of four dressed (anti-)quark propagators. The pole 
residue of the scattering equation is the homogeneous BSE 
for the four-quark Bethe-Salpeter (BS) amplitude shown in 
Eig. 1 . In compact notation it is written as 

r = KGo r, (2) 

which has a solution only if the T matrix has a pole. The 
exact kernel K is the sum of two-, three- and four-body 
irreducible interactions. In Eig. 1 we already use an ap¬ 
proximated kernel, where we omitted all irreducible three- 
and four-body interactions (the resulting equation can be 
rewritten as a Eaddeev-Yakubovsky equation [18]). 

This approximation is severe and can be disputed. Cer¬ 
tainly, in a more complete approach, these terms need to 
be analyzed in detail and taken into account accordingly. 
The main reason to omit these terms here is simplicity. As 
will become clear below, solving the relativistic four-body 
equation is a tremendous numerical task and some simpli¬ 
fications have to be made in order to make the equations 
tractable. However, we also believe that these approxima¬ 
tions receive an important justification a posteriori. As al¬ 
ready stated in the introduction, we And that the structure 
of the tetraquark wave functions are strongly dominated by 
the meson poles in the two-body interaction channels. As a 
consequence, the masses of the resulting bound states rely 
almost completely on the location of these poles. This fea¬ 
ture of the four-body BSE will not be changed by three- 
body or four-body interactions, unless these introduce new 
analytic structures into the equation that have a stronger 
influence on the integration region of the BSE than the 
meson poles. No such structures are known in the three- 
(anti)quark and four-(anti)quark channels. We therefore 
believe that the emerging ‘meson molecule’ picture, in the 
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sense explained in the introduction, is robust with respect 
to this approximation. 

In addition, we note that the introduction of irreducible 
four-quark interactions a priori does not have any prejudice 
towards a dq-dq or meson molecule picture. For exam¬ 
ple, the instanton-induced six-fermion effective interaction 
used in Ref. [14] introduces a tetraquark-gg mixing term 
that has (anti-)diquark channels for the tetraquark states but 
also terms where the tetraquark is in its meson-meson fla¬ 
vor state. Pending detailed investigations, this must be a 
general feature of the four-body kernel. 

Coming back to Fig. 1, the speciflc form of the remaining 
two-body interactions is necessary to prevent overcount¬ 
ing [16, 19, 20]. The equation in the figure is written in the 
dq-dq topology (12)(34); there are two further permuta¬ 
tions (23)(14) and (31)(24) with meson-meson topologies. 
The two-body kernels that appear in Fig. 1 must be consis¬ 
tent with the underlying quark-gluon structure to preserve 
QCD’s chiral symmetry. This is achieved by employing 
a rainbow-ladder kernel, where the qq and qq interaction 
is generated by iterated non-perturbative gluon exchange. 
The Dyson-Schwinger equation for the quark propagator is 
solved with this interaction. All relevant formulas and in¬ 
put values are given in Refs. [21, 22] where the approach is 
applied to baryons; we will not repeat them here for brevity. 
The rainbow-ladder setup describes the phenomenology of 
pseudoscalar and vector mesons as well as baryon octet 
and decuplet ground states reasonably well; see [23] for 
reviews. This implies not only mass spectra but also their 
form factors and other properties, and it extends to char- 
monium and bottomonium spectra [24]. 

One should note that rainbow-ladder alone does not pro¬ 
vide satisfactory results for scalar and axialvector mesons 
(the ‘p waves’ in the quark model). Viewed as a 
qq state, rainbow-ladder produces a d—meson mass of 
600 ... 700 MeV [25], with potentially sizeable correc¬ 
tions beyond rainbow-ladder [26] contrary to the masses of 
pseudoscalar and vector mesons. In any case, since these 
calculations support ideal flavor mixing they still lead to 
the mass ordering {d, Oo} —s- /t —)• /o, which means 
there must be further effects at play. On the other hand, 
a scalar qqqq state does form an s—wave orbital ground 
state. Hence, employing a rainbow-ladder kernel in the 
four-body equation (as opposed to a scalar qq meson) is 
well motivated from the aforementioned meson and baryon 
studies and its reliability can be judged from their results. 

Tetraquark amplitude. — The 0+ tetraquark BS am¬ 
plitude has the form Tap-y&ip-, q, k, P), where Greek sub¬ 
scripts are Dirac indices, P is the total momentum with 
P^ = —M^, and M is the tetraquark mass. A convenient 
choice for the three relative momenta are the Mandelstam 
momenta in the meson-meson and dq-dq configurations. 


i.e., the respective sums of the pair momenta: 

P =\{V2+PZ -Vl-Vi) , 4 

Q = l(P3 +Pl -P 2 -P 4 ) , ^ = ^Pi- (3) 

k =\{P1+P2-P3 -P 4 ) , 

Apart from color and flavor, the Dirac part of the full am¬ 
plitude is constructed from 256 tensor structures: 

256 

T{p, q, k,P) = Y^ /„(...) Tnip, q, k, P). (4) 

n—1 

The BSE leads to a system of coupled integral equations for 
the dressing functions /„ which depend on nine Lorentz 
invariants: the ‘Mandelstam variables’ q^, k^ and six 
further angular variables (p ■ q, p ■ P, etc.). This is also the 
main difficulty in solving the BSE numerically: even with 
only 10 points for each Lorentz invariant, each dressing 
function would be defined on 10® grid points. 

To make the problem tractable, we arrange the Lorentz 
invariants into multiplets of the permutation group S 4 [27]. 
This allows us to switch on groups of variables separately 
in the solution process and thereby judge their importance 
without destroying the symmetries of the system. The 
Mandelstam variables can be grouped into a symmetric sin¬ 
glet Sq and a doublet P, 

\/3(g® -p2) 

p^ + q'^- 2A:" J ’ 

(5) 

whereas the remaining six angular variables form two 
triplets 7i, 72- Hence, the dressing functions can be writ¬ 
ten as /„ = /ra(‘5o, 77 , 7i, 72)- The doublet phase space 
that remains invariant under the equation forms the interior 
of a triangle bound by the lines p® = 0, = 0 and = 0 

(see Eig. 2), whereas the triplets form a tetrahedron and a 
sphere, cf. Ref. [27]. 

To construct the Dirac tensor basis, we orthogonalize the 
momenta p,q,k, P to obtain four orthonormal momenta 
nf (z = 1... 4) that are mutually transverse. The set 

{1, (6) 

with i < j < k < I consists of 16 elements; commutators 
are not necessary because = —ifi. Taking all ten¬ 
sor products of Eq. (6) with itself yields 256 independent 
tensor structures. On the other hand, with the help of 

i, ik ii = -^75 , e := nj , (7) 

where {ijkl} is an even permutation of {1234}, all ele¬ 
ments with three or four slashes can be reduced to those 
with two at most, which leaves eight elements in Eq. (6). If 
we write Oi = 1 and O 2 = £ 75 . identify = P^ with 
the normalized total momentum, and express ^. in terms of 
the positive/negative-energy projectors A± = (1 ± j'l, 
then a complete, orthonormal, 256-dimensional positive- 
parity basis for the BS amplitude is given by 

r„(p, g, fe, P) = Fi Aa TsC 0 ( 7^75 F^ (8) 


^ _p^ + g^ + A:^ 1 

4 ’ ^- 4^0 
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FIG. 2. Top: Mandelstam triangle whose axes (not shown) are the 
doublet variables in Eq. (5). The top/right/left edge of the trian¬ 
gle corresponds to vanishing k^, or q^, respectively. The con¬ 
tour plot is representative and shows the magnitude of the leading 
dressing function. Bottom: Same plot in polar coordinates, where 
the radius 0 < r < 1 measures the distance from the center to the 
edges. The gap is due to missing diquark poles. 

with Fi G {1, We inserted the comhina- 

tion 75(7 <8) (7^75 for the dq-dq topology (12)(34); all fur¬ 
ther structures such as 7^(7 (8> (7^7'' hut also those in the 
meson-meson topologies are linearly dependent. 

In practice we are interested in partial-wave bases whose 
elements are eigenstates of the total quark spin and orbital 
angular momentum in the tetraquark’s rest frame. The con¬ 
struction is analogous to the partial-wave analysis of the 
nucleon’s Faddeev amplitude [22]. Since the tetraquark 
carries total angular momentum J = 0, the spin eigenvalue 
coincides with the orbital angular momentum (s = 1). The 
analysis yields 32 s waves, 144 p waves and 80 d waves. 
To obtain the s waves, replace F^ (8) F^ in Eq. (8) by 1 (8) 1 
and ^ 7 ^ ( 8 > Jp, where 7 ^ = 7 ^ ~ ^4 is transverse to 
the total momentum. This gives 2 x 16 = 32 tensor struc¬ 
tures. For the calculation we keep the 16 dominant s waves 
that do not depend on any relative momentum, i.e., those 
with G {1, 75}. This is a closed set under Fierz 

transformations: it contains dq-dq structures such as those 
formed by scalar and axialvector diquarks, 

i7^CU(C^7n,s, ... 

but by Fierz identities it automatically also includes all 


meson-meson structures in the crossed channels: 

( 75)07 1 (7 )q :7 (7 )/3i5 ) Stc. 

The associated error of neglecting the p— and d—wave ten¬ 
sors is discussed in the results section. 

Concerning color and flavor, the SU (3) decomposition 

3(8)3(8)3(8)3 = (3©6)(8)(3(8)6) 

_ (9) 

= 1 © 1 © 8 © 8 © 8 © 8 © 10 © 10©27 

provides two singlets via 3 © 3 and 6 © 6 (or, in the 
meson-meson configurations, 1 © 1 and 8 © 8 ). The com¬ 
plete list can be found in [28, 29] or reconstructed from 
the Clebsch-Gordan tables in Refs. [30, 31]; alternatively, 
the program of Ref. [32] can be used. Even though the in¬ 
teraction in the color-sextet diquark channel is repulsive, 
we take into account both color singlets because they are 
necessary for Fierz completeness. Regarding flavor, all of 
the above multiplets may contribute in general and those 
with same isospin and hypercharge can mix. In our case 
we can simply factor out the flavor structure because the 
rainbow-ladder kernel is flavor-blind; different flavor states 
will show up as excited states in the BSE solution. 

Results. — In practice the BSE is solved as an 
eigenvalue equation, schematically written as KGoTi = 
Xi{P^)Ti. The largest eigenvalue describes the ground 
state and the remaining ones the excited states. Upon cal¬ 
culating the eigenvalue spectrum, \i{P^ = —M^) = 1 
recovers the original equation and thereby allows one to 
extract the mass and BS amplitude of the state. 

The multiplet analysis around Eq. (5) provides us with a 
convenient tool to isolate the relevant momentum regions 
probed by the equation. The largest eigenvalue for the sys¬ 
tem consisting of four light quarks is shown in Fig. 3. The 
second curve from the bottom shows the result obtained 
by restricting the momentum dependence to Sq only, i.e., 
/„(5o,F>,7I,7^) /„(5o). The corresponding bound 

state mass is M ~ 1.5 GeV, which is essentially ‘four 
times the constituent-quark mass’. Taking into account ei¬ 
ther of the triplets does not change this behavior signifi¬ 
cantly. However, implementing the doublet has a drastic 
effect: the eigenvalue is now almost flat and close to 1, 
with a crossing at ~ 350 MeV. 

Therefore, it is the doublet phase space that carries the 
relevant momentum dependence. This can be understood 
from Fig. 1: upon iteration, the BSE will generate two- 
body qq and qq scattering matrices with their own (meson 
and diquark) singularity structures. These two-body poles 
appear at timelike values of the Mandelstam variables, i.e., 
in the exterior of the doublet triangle in Fig. 2. The onset of 
the pion poles at <0 and q^ < 0 can be clearly seen in 
the contour plot as well as in the 3D plot of Fig. 2: the pion 
mass is small and the pion poles are close to the triangle, 
whereas the diquarks have no visible effect because their 
mass scales are much larger: nisc ~ 800 MeV and ~ 
1 GeV in rainbow-ladder [33]. 
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FIG. 3. Largest BSE eigenvalue for the tetraquark with four light 
quarks after switching on different multiplets. The band in the 
inset shows the result of linear and quadratic extrapolations; the 
resulting vertical bar gives the extracted a mass. 


This feature is what drives the tetraquark mass towards 
low values: although the four-quark equation knows a pri¬ 
ori nothing about mesons or diquarks, it generates their 
pole structures dynamically and these influence the hehav- 
ior of the dressing functions from the exterior of the inte¬ 
gration domain. Since the system is driven hy the would-he 
pseudoscalar Goldstone hosons, the tetraquark mass can he 
understood as a remnant of spontaneous chiral symmetry 
breaking. On the other hand, this also naturally leads to 
a meson-molecule interpretation (in the sense explained in 
the introduction) simply because the diquark mass scales 
are too large to be relevant. 

The smallness of the pion mass also generates a physi¬ 
cal threshold: if M > 2m^, the poles enter the integration 
domain and the tetraquark becomes a resonance. In that 
case the tetraquark pole in the four-quark scattering matrix 
will disappear from the real axis and ultimately one should 
solve the BSE in the complex plane, also taking into ac¬ 
count the poles in the integration domain. This is, however, 
a rather formidable task which has not even been accom¬ 
plished in simpler systems so far. Hence, for the moment 
we resort to simple extrapolations to estimate the real part 
of the mass, although we are aware that crossing thresholds 
leads to non-analyticities [34] and the results of extrapola¬ 
tions should be interpreted with great care. 

To minimize numerical artifacts due to the arising pole 
structures, we approximate the dressing functions of the 
tetraquark BS amplitude by 




3 

i: 

3 




^ (gj + - {q, ■ py ’ 


( 10 ) 


where qj G {p, q, k} are the three Mandelstam momenta, 
rUj are the corresponding meson and diquark pole masses 
that appear in the doublet, and R is the radius of the tetra¬ 
hedron which forms the triplet Ti- The denominator is the 


M [GeV] 


strange 


charm 



FIG. 4. Left panel: tetraquark masses at fixed strange quark mass 
and varying light quark mass. The bands contain the extrapola¬ 
tion errors and the dotted curve is the two-pion threshold. The 
vertical dotted line shows the physical u/d mass and the vertical 
bars are the experimental values for the masses [4]. Right panel: 
quark-mass dependence of the tetraquark mass with four identical 
quarks. The band is the present result from the four-body equa¬ 
tion with extrapolation errors, compared to that obtained with 
coupled meson-meson/dg-dg equations (dash-dotted curve) [16]. 
The dotted curves are the respective thresholds. 


product of two-body poles and provides an accurate repre¬ 
sentation of the momentum dependence in the Mandelstam 
plane (Fig. 2); the pole residues are determined dynami¬ 
cally in the BSE solution. In addition to the Sq dependence, 
the numerator also captures the subleading effects provided 
by the tetrahedron that are visible in Fig. 3. 

The resulting current-quark mass dependence of the a 
mass is shown in the right panel of Fig. 4. The tetraquark 
is a resonance just slightly above tttt threshold; it only 
becomes a bound state in the charm-quark region. The 
mass is also larger than our earlier result in Ref. [16] 
obtained with the coupled meson-meson/dg'-dg equations 
(also shown in Fig. 4). However, this might be due to our 
omission of the remaining p— and d—wave tensor struc¬ 
tures which are likely to provide ~ 10% effects in analogy 
to three-quark systems [21]. Analogous conclusions apply 
for the multiplet partners of the a, shown in the left panel 
of Fig. 4; the k is dominated by iTvr poles and ao//o by 
ptr and KK poles. Our results at physical current-quark 
masses are 


= 0.348(13) GeV (11) 

together with 

M, = 0.750(12) GeV, 

= 1.081(28) GeV, ^ 

where the errors are due to the chosen extrapolation proce¬ 
dure (linear or quadratic). 
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Summary and conclusions. — To summarize, we 
presented a numerical solution of the four-quark Bethe- 
Salpeter equation for scalar tetraquarks, implementing a 
well-established rainbow-ladder (gluon-exchange) interac¬ 
tion. The resulting masses for the a, k and a ^/are 
in the ballpark of their experimental values. Barring fur¬ 
ther mixing effects, this provides support for a tetraquark 
interpretation of the light scalar mesons. Our calculation 
also illustrates how a dynamical generation of resonances 
emerges from the quark level: pseudoscalar-meson poles 
are generated in the solution process and dominate the dy¬ 
namics of the system; if they enter the integration domain 
the tetraquark becomes a resonance. In that sense these 
tetraquarks are not diquark-antidiquark states but predom¬ 
inantly ‘meson molecules’, which explains their mass or¬ 
dering and their decay widths into the hadronic channels 
as an indirect consequence of QCD’s spontaneous chiral 
symmetry breaking. 

While we believe that our calculation makes a tetraquark 
interpretation of the light scalar meson sector further plau¬ 
sible, much can be done to further test and improve the 
framework. On systematic grounds it would be highly in¬ 
teresting to study tetraquarks with other quantum numbers. 
Within the currently employed framework and its approx¬ 
imations this is feasible in principle, but it requires huge 
computational resources. Furthermore, one would like to 
lift some of the approximations made so far. Expanding 
the framework to include also the p— and d—wave ten¬ 
sor structures is conceptually straightforward, but, again, 
cost-intensive in terms of CPU time. The introduction of 
additional three-body and four-body irreducible forces re¬ 
quires some conceptional effort; work in this direction is 
underway. 

An important future application of our approach are 
tetraquarks in the charmonium region. One of the most 
interesting problems there is to identify the internal struc¬ 
ture of potential candidates, i.e. to decipher whether 
they are dominated by diquark-antidiquark (nc-nc), meson 
molecule (nc-nc) or hadrocharmonium (nn-cc) configura¬ 
tions. These correspond to poles along the edges of the 
Mandelstam triangle in Fig. 2. Our approach has the po¬ 
tential to answer this question dynamically from the under¬ 
lying quark-gluon interaction, thus bypassing the assump¬ 
tions made by many models. 
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